Charged State of a Spherical Plasma in Vacuum 
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The stationary state of a spherically symmetric plasma configuration is investigated in the limit 
of immobile ions and weak collisions. Configurations with small radii are positively charged as a 
significant fraction of the electron population evaporates during the equilibration process, leaving 
behind an electron distribution function with an energy cutoff. Such charged plasma configurations 
are of interest for the study of Coulomb explosions and ion acceleration from small clusters irradiated 
by ultraintense laser pulses and for the investigation of ion bunches propagation in a plasma. 

PACS numbers: 



I. INTRODUCTION 



The interaction of ultraintense laser pulses with solid 
targets leads to the formation of plasmas with unusual 
properties in terms of particle energy distributions and 
of spatial mass and charge density distributions. Such 
properties can be exploited in order to obtain sources 
of high energy electromagnetic radiation and charged 
particle beams with unprecedented intensities and time 
and space resolutions. For the intensities of present day 
ultrashort, superintense laser pulses, the energy that the 
ions in a target acquire due to direct interaction with the 
electromagnetic fields of the laser pulse is usually small, 
while the energy of the plasma electrons can be of the 
order of tens of KcV. These "hot" electrons expand until 
their "pressure" is balanced by the electrostatic field 
that sets up due to spatial charge separation Q, |E Q- 
This process leads to a steady state configuration (SSC) 
which is reached in a time of the order of some electron 
plasma periods. Thus, since the ion response time is 
much longer than that of the electrons, SSC can be 
achieved before the ions can depart significantly from 
their initial configuration. Afterwards ion acceleration 
takes place, as predicted theoretically 0, IE IE EL IE IE IE 
IIJJlEJSjQjEJii^ and crafimed experimentally 

EE EE E3, EE Ha Hil Ha m, in, HE Hal. cieariy, 

this description docs not apply to highly relativistic 



regimes like those described in 
The aim of this paper is to present a combined analytical 
and numerical investigation of three dimensional SSC 
characterized by a hot electron plasma and a cold 
(immobile) ion core. Such configurations are especially 
appealing because, contrary to one dimensional ge- 
ometries, in three dimensional cases charged SSC are 
expected to set up. In fact, while in the former case an 
infinite energy is required in order to bring a charge to 
infinity, in the latter the energy necessary for electron 
evaporation is finite. In particular, we show that the 
SSC charging up and the energy distribution of the 
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electrons depend on the ratio between the radius of 
the ion core and the electron Debye length and on the 
history of the electron expansion. The understanding 
of the SSC charging up with immobile ions and of 
the electron energy distribution is relevant to many 
experimental conditions as these processes affect the 
way in which ions are accelerated on longer time scales 
when the constraint of a fixed ion core is removed. In 
the case of the Coulomb explosion of a small cluster 
EE HE EE EE IH EE these processes affect the 
value of the maximum energy that the ions can acquire 
in the acceleration process. Furthermore, in applications 
related to proton imaging [3^ . l4fl l4lj and to the prop- 
agation of ion beams in solid targets EE|j the Coulomb 
repulsion and the screening effect of the electrons can 
strongly affect the proton trajectories. This is also 
the case for applications of proton laser acceleration to 
hadron therapy HE E3, EE llE E3, EE EE EE| , where a 
very precise collimation of the proton beam and a high 
energy resolution are essential. Moreover, as discussed 
in |52| i the topic of the charging of a spherical plasma in 
less extreme conditions can play a key role in the study 
of dusty plasmas |sE| > an d aerosols |5J| . 
The paper is organized as follows. In Sec.^ simple 
relationships are derived analytically on the basis of 
two schematic models that are introduced in order to 
grasp the main features of the plasma charging process 
in different collisionality regimes. In Sec. lIIII we present 
the results of a series of numerical simulations obtained 
with a one-dimensional Particle in cell code (PIC) in a 
spherical geometry and with fixed ions. We then com- 
pare the numerically obtained charge values and electron 
energy distributions with those obtained from the analyt- 
ical models. Finally, conclusions are presented in Sec . 11 VI 



II. SIMPLIFIED MODELS OF THE CHARGING 
PROCESS 

In this section we discuss two simplified models that al- 
low us to identify the main physical features of the plasma 
charging process. These models rely on assigning a sim- 
plified condition for the electrons to leave the ion core 



2 



and on two different rules for the energy redistribution 
of the remaining electrons. 

As a starting configuration we assume the following elec- 
tron and ion density profiles 



n e (r) = ni(r) = n 9(1 - r/R), 



(1) 



with 6(x) = 1 for x > and 8(x) — for x < 0. Here R is 
the radius of the ion plasma core and Nq = no(47ri? 3 /3) 
the ion and electron initial particle number. We denote 
by N e the time dependent number of electrons inside 
the ion core. 

Initially electrons have a Maxwellian energy distribution 
with temperature To. Moreover, in these models the 
electron density is taken to be uniform inside the ion 
core. As a further simplification, we assume that, on 
average, the radial crossing of the electron trajectories 
does not lead to a relative redistribution of the charge 
in front and behind each electron outside the core. This 
simplification allows us to assume that, after leaving the 
ion core each electron moves as if in an effective time 
independent Coulomb field. Hence, the condition for an 
electron to reach infinity is that it has a positive total 
energy when it reaches the ion core surface at r = R. In 
this model the particles which satisfy this condition are 
said to "evaporate" and are assumed to be lost when 
they reach r = R. On the contrary, the electrons that 
have a negative energy at r = R are assumed to remain 
inside the ion core. Furthermore, we assume that inside 
the ion core the electrons move as free particles. The 
evaporation of the electrons with positive total energy 
at r = R changes the number of electrons N e , the total 
energy of the system inside the ion core and causes an 
energy redistribution of the remaining electrons. 
We discuss the "collisional" regime and the "collision- 
less" one. In the first one, the electrons which are not 
evaporated thcrmalize at a temperature T, which turns 
out to be a decreasing function of time. In the second 
regime no thcrmalization occurs, and the evaporation 
causes a progressive depiction in the high energy tail 
of the electron distribution function, which remains 
isotropic in velocity space. 

In what follows, lengths will be measured in units of the 
initial Debye length \d = (T /47rnoe 2 ) 1 / 2 , with e the 
absolute value of the electronic charge, time t in units 
of uj~e = (47re 2 no/TO e ) -1 / 2 , with m e the electron mass, 
energies in units of the initial electron temperature To, 
velocities in units of the initial electron thermal speed 
Vthx) = y/To/m e , mass in units of the electron mass and 
particle numbers in units of Nq. Since inside the ion 
core the electron density is taken to be uniform, with the 
adopted normalization the normalized electron density 
n e and the normalized total number of electrons N e are 
numerically equal. 



A. Collisional regime 

If the electrons inside the ion core are re-thermalized 
by collisions, their velocity distribution function remains 
Maxwellian and their time dependent kinetic energy is 
given by Uk = 3N e T/2. The electron evaporation rate is 
obtained by calculating the flux of electrons with positive 
total energy through the core surface. Defining the time 
dependent quantity <p T = e$ R /T = (l-N e )R 2 /3T, with 
&r the electrostatic potential at the ion core surface, the 
positive total energy condition corresponds to v 2 /2T > 
d>T- Thus we obtain 



dt 



0- + <h) 



/2tt 



> T N 



(2) 



where r = R/\T is the electron crossing time inside 
the ion core. Analogously, the energy flux $u of the 
evaporating particle can be computed by noting that 
each evaporating electron carries away the residual en- 
ergy v 2 /2T — <Pt- Then we obtain 



4>, 



3 (2 + T ) 



N„T. 



(3) 



The total energy of the system can be written as U = 
3N e T/2 + £/$ where U$ is the electrostatic energy of 
the charged configuration which increases as the electrons 
evaporate as 



(4) 



Thus, from the total energy balance we obtain for the 
time change of the kinetic energy of the system 



d(3N e T/2) 



3 {2 + 4> T ) 



dt 



e -<p T N T 



(5) 



which, coupled to Eq. J2J), gives the time evolution of 
the temperature T. 



B. Collisionless regime 

If on the contrary we assume that plasma elec- 
trons inside the ion core are not significantly affected 
by collisions, their distribution function becomes non- 
Maxwcllian. We assume that the electron distribution 
remains homogeneous in coordinate space and isotropic 
in velocity space. Thus, denoting by Ng the time de- 
pendent number of electrons with kinetic energy (nor- 
malized on the initial tempcratute T ) in the interval 
[£,£ + dS], and introducing the time dependent quantity 
4>o = e^ii = |(1 — N e )R 2 , which differs from 4>t in the 
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previous by the normalized temperature factor 1/T. we 
obtain 



dN £ 
dt 



3 VE 

2V2 R 



0(8 -fa) N £ . 



(0) 



This implies that the evaporation of the electron popu- 
lation with energy £ stops at a well defined time t = t £ , 
where t £ is such that £ = (f>o(t £ ). Therefore, 



N £ = N £ (0)e- t/td for t<t £ 
= N £ (0)e- t£/td for t>t £ 



(7) 



with td = the ^-dependent decay time 

and N £ (0) the electron kinetic energy distribution at the 
initial time t = 0. We assume the electron velocities at 
time t = to be Maxwellian distributed, hence the initial 
electron kinetic energy distribution N £ (0) is given by 



N £ (0)=N £ , o = ^e~ £ V8. 

\/7T 



(8) 



The electron number N e , and therefore <po, can thus be 
calculated performing numerically, at fixed t, the inte- 
gral N e = J d£N £ . Note that in this collisionlcss model a 
rough estimate of the asymptotic electron number could 
be obtained by approximating the final electron distri- 
bution function with the initial one for £ < 8*, and 
with zero for £ > £*. We can then determine the cut- 
off energy 8* self consistently by equating its value to 
the electrostatic energy of the configuration with charge 
Q{8*) = J™N £ (0)d8, 



£* = [Q{£*) i? 2 /3] 



C. Discussion of the analytical models 



(9) 



Numerical integration of Eqs. J2J[SJ shows that both in 
the collisionlcss and in the collisional regime the electron 
evaporation rate is significantly reduced from its initial 
value when the quantities </>t,</>o become of order unity. 
This means that, in the collisional regime, the evapora- 
tion nearly stops when the electron electrostatic energy, 
which is an increasing function of time, is of the order 
of the electron temperature which decreases with time. 
Whereas in the collisionless case the electron evaporation 
is significant only until the electron electrostatic energy 
is of the order of the initial average electron kinetic en- 
ergy. 

The initial evaporation rate, which is obviously the same 
in both, the collisional and in the collisionlcss case, is 
obtained from either Eq. J5J) or Eq. 10 which give 



dNe 

dt 



t=o 



1 



/2tt R 



(10) 



A linear approximation of N e is feasible as long as 
i<f,t where 



is the initial characteristic evaporation time and i is de- 
fined in the two different collisionality regimes by cither 
the condition 4>t\i ~ 1 or 4>o\i ~ 1 an d which can be 
roughly evaluated as 



t = 3t*/R 2 



(12) 



2ttR/3 



(11) 



Therefore, if t* < t, i.e., for ion core radii satisfying 
R < V3, the charging process continues until almost 
complete depletion of the electron population. 
Finally, we note that the time dependent electron energy 
distribution predicted in the collisionless regime is highly 
non thermal. The fact that the electron evaporation 
only occurs for those particles with kinetic energy £ 
satisfying £ > </>o| t causes a depletion of the high energy 
tail of the electron distribution function, as will be 
examined in detail in Sec. IIIII 



III. PIC SIMULATIONS AND COMPARISON 
WITH THE ANALYTICAL RESULTS 

Our PIC simulations are initialized with the electron 
and ion density profiles n e ,rii given by Eq. The 
initial electron distribution function is Maxwellian with 
temperature To. At time t — the electrons are allowed 
to move. During their expansion, the electrons that reach 
the border of the simulation box with positive total en- 
ergy are removed. Therefore, as the total number of elec- 
trons decreases with time, the plasma acquires a net posi- 
tive charge and an electrostatic potential sets up. Finally 
a stationary state is reached where no more electrons 
reach the boundary. We denote by N 00 ,n 00 ,E 00 ,N £t00 
the electron number, the electron density, the electric 
field profile, and the electron kinetic energy distribution 
once SSC has been reached. As expected, our simula- 
tions show that these quantities depend on the ion core 
radius R. 

The results regarding the profile of both the electric field 
Eoc and of the electron density Hoc for two different ion 
plasma radii R are presented in Fig.^ As one can see 
in Fig.^D, the electrons which are outside the ion sphere 
are confined in a region whose width is of the order of a 
few Debye lengths. 

The numerical results regarding the electron kinetic 
energy distribution in the stationary state are presented 
in Fig.|2 In the figure, the ratio N £)00 /N £i o, is shown 
versus the electron kinetic energy £ in semi-logarithmic 
scale, for several values of the ion plasma radius R. These 
results show that the electron kinetic energy distribution 
is highly non-thermal. It exhibits a cut in its high energy 
tail, and the cutoff energy increases with the ion plasma 
radius R. 

Regarding the comparison with the semi-infinite, pla- 
nar case, our results show that both the electric field 
and the electron density are very similar to those pre- 
sented in |(| only as long as the ion core radius is greater 
than several tens of Debye lengths. Since, contrary to 
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FIG. 1: Spatial profile of the electric field Eoo (panel (a)) 
and electron density n x (panel (b)) for R = 10 (dashed line) 
and R = 40 (solid line). For comparison they are plotted 
together with the profile that would be obtained in the semi- 
infinite case (dotted lines) 0. The straight line in panel (a) 
corresponds to the maximum value of the electric field in the 
semi-infinite case, i.e, E = yj2/e. 
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FIG. 2: Behaviour of the ratio between the final and the initial 
electron kinetic energy distributions for different ion plasma 
radii. 



FIG. 3: Comparison between the stationary state value of 
the dimensionless charge inside the ion core predicted by the 
collisional model (triangles), the collisionless model (squares), 
the PIC simulations (stars) , and the numerical solution of Eq. 
© (solid line). A magnification of the region 15 < R < 45 is 
also shown in the figure. 



a onc-dimcnsional configuration, in the case of a three- 
dimensional configuration the energy required for the 
electrons to evaporate is finite, the differences observed 
are mainly due to the cutoff in the electron high energy 
tail. Such cutoff is responsible for the electron density 
depletion observed outside the ion core (see Fig.^j) and, 
consequently, for the corresponding electric field profile. 
As shown in Fig.[TJi, in the limit R 3> 1 the value of the 
dimensionless electric field at the surface of the ion core 
is almost independent of R, thus the net dimensionless 
charge confined within the ion core scales approximately 
as 1/R. 

The comparison between the charge value Qoo = 1 — Noo 
obtained numerically and that predicted analytically is 
shown in Fig. [3 for different values of the ion core ra- 
dius R. Note that the value of Qoo obtained numerically 
includes all SSC electrons, i.e., the electrons inside the 
ionic sphere and those in the surrounding halo. It is seen 
that the agreement among the numerical results and the 
values obtained in the collisionless regime is very good 
in the whole range 5 — 40. With regard to the thermal 
model adopted for the collisional regime we remark that 
for small radii it predicts a moderately larger value of 
Qoo, but the two different regimes lead to very close val- 
ues Qoo in case of large radii. 

Our results indicate that collisions can affect the charg- 
ing up process only for small ion core radii. This result 
can be explained by noting that the potential due to the 
electron expansion scales as R . Thus, as the potential 
barrier increases, the fraction of the electron population 
that, because of Coulomb collisions, reaches a positive 
total energy and can thus leave the ion core decreases. 
However, the collisional thcrmalization of the distribu- 
tion function assumed in Sec . III Al can only be expected 
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FIG. 4: Comparison between the cutoff energy £ cut as a func- 
tion of the ion core radius R predicted by the collisionless 
model (squares) and by the PIC simulations (stars). 



to to apply when the ion core radius is much larger than 
the electron mean free path, whereas in most plasma con- 
ditions the mean free path due to Coulomb collisions is 
much greater than the Debye length. 
It is worth to note that a precise fit of the numerically- 
obtained plasma charge state is given by the following 
Pade approximation 



1 + aR 



1 + bR + cR 2 ' 



(13) 



with a = 5.6 ■ 1CT 3 , b = 1.4 • 1CT 1 and c = 5.5 ■ lO" 2 . 
With regard to the high tail of the electron kinetic energy 
distribution we define the cutoff energy £ cut as the en- 
ergy satisfying the relation (N e ,oo/Ne,o)\e=e cut = 10 _1 . 
The analytical and numerical results obtained for £ cu t 
are shown in Fig. 21 for different values of the ion core 
radius R. Taking into account these results, we remark 
that the collisionless model, although not capable of re- 
producing the fine details of the whole electron kinetic 



energy distribution observed in the simulations, never- 
theless it predicts the cut in the distribution high energy 
tail with great accuracy. 



IV. CONCLUSIONS 

In this paper we have investigated the charging up pro- 
cess of a spherically symmetric plasma configuration in 
vacuum in the limit of immobile ions. Two different sim- 
plified models have been presented. With the help of 
these models we have established scaling laws relating 
the steady state total charge and electron energy distri- 
bution on the radius of the ion core, normalized in terms 
of the initial electron Debye length. These scalings have 
been validated by mean of spherical one-dimensional PIC 
simulations. The agreement we find is overall very good. 
Charged plasma configurations such as those investigated 
in this paper occur naturally in the interaction of ultrain- 
tense laser pulses with matter and are of interest, e.g., 
for setting the initial conditions in the study of Coulomb 
explosions and ion acceleration from small clusters irra- 
diated by ultraintense laser pulses. 

In particular, regarding the problem of cluster expansion 
the analytical and numerical results that we have pre- 
sented show that a spherically symmetric configuration 
of cold ions and hot electrons, which is the typical start- 
ing configuration in cluster expansion experiments, does 
not evolve towards a neutral configuration, in contrast 
with the onc-dimcnsional planar case. This charging up 
effect can strongly modify the maximum energy that the 
ions can gain and the typical timescale on which their 
acceleration occurs. 
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